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Abstract. In this paper we show the existence of a large class of spherically 
symmetric data d (on a spacelike hypersurface S), from which a perfect fluid 
spacetime (surrounded by vacuum) develops. This spacetime contains an event 
horizon (with trapped surfaces behind it). The data d are regular and innociuous, 
i.e. the data— surface S does not contain any point of the horizon or of the trapped 
surface area. The occurence of the horizon (and trapped surfaces) is stable under 
small spherically symmetric variations of the data d. 

We give auxiliary data on an auxiliary hypersurface H and also on the star 
boundary; then we solve Einstein's equations for perfect fluid in the future and past 
of H . Our solution induces the above mentioned data d on some chosen spacelike 
hypersurface S in the past of H . By construction H turns out to be the matter part 
of the horizon, once we attach a vacuum to our matter spacetime. Obviously, from 
these data d on S it develops (into the future) the event horizon H . 

We solve the constraint equations for the auxiliary data posed on the null- 
surface H. This reduces the choice of these data to the choice of the density p 
and R := [the orbit curvature] Our data fulfil positivity of p, 2m /R = 1 (at the 
star boundary) and other properties. This is archieved by an algorithm, which for 
given p yields R (from an input parameter function h G ([0, 1] , ]0, — oo[)). 
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1. Introduction 

Let us confine ourselves to spherically symmetric spacetimes. In many interesting 
papers (e.g. [1],[2]) there was shown that sufficient mass-concentration induces the 
forming of trapped surfaces: certain conditions have been posed on data d given on 
some spacelike hypersurface S; it was shown that these conditions imply the existence 
of trapped surfaces on S. 

One would like to extend these results into three directions 

(i) The data on S should be mnoctuous, i.e. the trapped surface should not occur 
on S, but rather in the future of S; 

(ii) One would like to show the forming of an event horizon (a global object) along 
with trapped surfaces (a local object); (see [3], where null-surfaces are used 
instead of the spacelike surface S) 
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(iii) One would like to show the existence of a (unique) spacetime, which develops 
from the data d. In [2] (and also in [3], [4]), however the existence of such a 
spacetime has been assumed. 

In the present paper we pose auxiliary data rffj on an auxiliary hypersurfacej H , 
(which will turn out to be part of the event horizon); furthermore, we pose data on the 
star boundary. Then we show the local (in time) existence of a (unique) solution]: Vi 
assuming the above data; Vi is a perfect fluid spacetime containing also a part, which 
lies in the past of H . Now, in the past of H we choose a spacelike hypersurface S; our 
solution Vi induces data d on S. From these data d on S evolves - into the future of 
S - a unique solution; this solution obviously is our above solution Vi, which contains 
the hypersurface H in the future of S. The null-extension of H into the vaccuum, 
which we attach uniquely to Vi, is (by construction) an event horizon. 

The above construction means the following: we have constructed data d (on 
some spacelike hypersurface S), from which an event horizon (with trapped surfaces 
in its future) evolves. The occurence of the event horizon is stable (under spherically 
symmetric variations of the data d). The properties (i),(ii),(iii) are fulfilled. However, 
the mass-concentration will not be controlled on S but on iJ , on which the auxiliary 
data (ifj are given. The mass-concentration is controlled as follows (see Remark 13 c), 
d)): 

^<lonF\5, ^ = lonFn5, (1) 

K K 

where B is the star boundary and m denotes Misner's mass§ scalar field (see [5]-[7]) 

— 1/2 

and R{P) '.= [curvature of the orbit (of spherical symmetry) at the point P] 
The condition (1) is sharp: if we would replace (1) by ^ > 1 on iJ , then H would 
not be a horizon (see also Remark 13 c), d)). 

It may be possible to extend the results of [2] such that (i) is fulfilled; Moreover, 
(iii) may also be satisfied in a neighbourhood of S (using existence theorems for local 
solutions). However, in order to fulfil the global property (ii), one would need some 
global\\ existence theorem on solutions of Einstein's equations for perfect fluid. It 
will be difficult to obtain such global existence theorem, since the fiuid can evolve 
shell-crossings or shocks. 

One can avoid such difficulties, if one couples Einstein's equations not to a perfect 
fiuid but to a scalar field, which has been done in [8]. A global existence theorem was 
proved in [8], and the properties (i),(ii),(iii) are also fulfilled. 

The first author is writing his doctorial thesis, which includes the material of this 
paper as well as some numerical investigations. For details we refer to this thesis. 

2. General assumptions, definitions, and coordinates 

We shall assume that our spacetime is spherically symmetric (with one timelike 
centreline) and that there is given a smooth equation of state p(p) with density p 
and pressure p, i.e. 

t see Remark 9 a). 

I Actually, we shall construct this solution Vi in two steps; (see Figure 1, where Vi is in fact Vp U Vh)- 
§ This mass has been introduced by Tolman [9]. 

II In this paper we avoid using global solutions by choosing our (auxiliary) data— surface i7 to be a 
rau//— surface. 



(indices fi,iy . . . run from 1 to 4) fulfilling 
dp 

< — < 1, j5(p) > for p > pb and p{pb) = 0, (2) 
dp 

where pb > is some constant! . 

We consider Emsiem's equations for perfect fluid 
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where denotes the fluid vector and g^" the metric field. We define Dt := 
u^dxf (dxi' '■= gfTr) and the unit vector field Df as being orthogonal to Df and 
to the orbits (of symmetry). Later, we shall consider a given (r, t)-coordinate plane 
((r,t) of Definition 2). For r = we define Df as its continuous extension, provided 
we deal with spacetimes endowed with a regular centreline at r = (e.g. in Lemma 
5). Moreover, we define the scalar fields 

R := l/\/[curvature of the orbits], U := DtR and F := DfR (3) 

Furthermore, we define the following auxiliary functions 

Ax 



<j){p):=-( ^Pj-^^ [{p(x) + x)] ^dx, n (/!>):= rib exp 



{p{x) + x) 



(4) 



where pb > (cf (2)), rib > are some constants. 

Lemma 1. a) Given and i?H as real-valued functions of r G [0, rb]. We define 

Bn{r):= e'-^^^^^^^^da;, An{r) := A% h{pn{x)) ■ R^^ix) Ax . (5) 

Jo Jo 

We assume that (ds'^ jU^dxi^ , p) is a spherically symmetric perfect fiuid spacetime. 
Then there exists a coordinate system such that 

ds^ = e^dr^ + i?2 (d«?2 + (sin^ ^) d^^) _ ^24, _^ 5/^ (^^) J^^)2 ^^q^ 
with A, R, (f> depending on {r,t) (a dash denotes ^) and with 

u^dxi^ = e~'^(9( , (i.e. co-movmg coordinates) (7) 

for r G [0, rb], t G [—T, T], where r = (and r = rb, respectively) denotes the centreline 
(and the star boundary, respectively). It follows that the scalar field 

A{r,t) := I A% h {p {x ,t)) ■ R^ {x ,t) ■ e^^'^'' ''^Mx depends on r only. (8) 



b) In addition we assume 

p = PH and i? = i?H on H := {(r,t) \t = 0,re [0,rb]}. (9) 
Now, we can (and shall) fix the coordinates of (6,7) by 

ct,(ry,,t) = ^ (PH (rb)) (for t G [-T, T]) A(r, 0) = An(r) (for r G [0, rj) (10) 

t We use the more convenient notations p{p) and ^ instead of p o p and p' o p. 
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The latter is equivalent to A = on iJ. It follows that iJ is a null-surface. 
Furthermore, (f) and A of (6) are calculated as 

4>ir,t) = 4>ipir,t)), e^/' =T-'{dr-Bi,dt)R. (11) 

Moreover (due to (8,10)), it holds 

A(r,t) = An(r). (12) 

Definition 2. We denote the coordinates fixed by (6,7,10) as {r,t)- coordinates. 

Remark on Lemma 1. Our (r, t)-coordinates are closely related to Misner's coordinates 
(see (13,14)). — Our (r, t)-coordinates are co-moving and spherically symmetric. It 
holds t = Q on H , which is convenient, when we shall give (/9h,-Rh) as ^data' on the 
surface H. These 'data' (pn, Rn) have been used for fixing the (r, t)-coordinates (see 
(10)). Incidentally, any quantity carrying an index H can be calculated from these 
functions (/9h,-Rh), see e.g. (43). These quantities depend on the variable r G [0,rb] 
only. — Note, that {(r,t) \ t = to} is in general not a null-surface, provided to ^ 0. 

Proof of Lemma 1. Let us denote by {r,t) the following coordinates used in [5]-[7]: 

ds2 = g,, (r, t) dr2 + ij2 (d^2 ^ (^^^^2 j^2) ^ ^_ ^ 

fixed byt A{r,t) = An{r), gjj(r^J) = -e^kPH(f,)) (^4) 

for r G [0,rb] with rb = rb and (t — Bii(r)) G [— T, T]. Our (r, t)-coordinates can be 
transformed into (r, t)-coordinates by 

t = t + Bn(r), r = r (5f = 5*, 5,- = 5, - 5(j5t) (15) 

Moreover, our scalar fields pi, R,U,T, A are the same ones as in [5]-[7] (same 
notation but different coordinate system). We obtain 

U = DtR = e-'*'dtR, V = DfR = e'^l'^ {dr - B'^dt) R (16) 
from (3), (6). One gains (11) from the analogous formula§ of [5]-[7] by (15). 

3. Einstein's equations in (r,t) coordinates 

Given and i?H as functions of r G [0,rb]. Using 4>,h, B]i, A]i of (4,5) we write 
Einstein's equations in (r, t)-coordinates as follows for r ^ 0: 



I) {r,t)- evolution equations, which are quasi-linear 
dtp = ■ (p + p) 



a' • 5{je'^(^) -T-^-drp-aA'rr-drU 
dp 



(17) 



which is called e in [5]— [7] 

Similiar results as in [5]— [7] have been obtained by [10]. 



dtU = a-^ 



• r 



2 '^.d^p_a^.B'^e^^P).f.T-drU 
dp dp 



dp R 



2R 



AirpR 
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(18) 
(19) 



dp 



dp WT 



" dp \ 2R 



where a := Ai:n{p)R^^ and a := F- 
abbreviations. 
II) {r,t)- con sir amis 

^_ 47rn(p)i?2 



dp 
AirpR 

2 



dp 
dp 



(20) 

Q-'Pip) are technical 



A' 



r2 _ 



2m 



m(r, t) 



47r [pS^i? - Bi^eMP^ ■ (p + p) ■ u] (x, t) ■ R^ {x, t) d. 



(21) 
(22) 
(23) 



The scalar field m is identical with m of Misner (cf [5]-[7]), provided the 
assumptions of Lemma 3 or Lemma 1 are fulfilled. Misner's m (r, t) denotes the mass 
(total energy) contained in the ball K (r,t) := {(x,t) \ x G [0,r]}, where he used the 
(r, t)-coordinates defined in (13), (14). 

Einstein's equations in (r, f)-coordinates are (Al) together with F = '^'^"^p^^ -dfR 

and F^ — [/^ — 1 + ^ = (see [6]). We transformed these equations via (15) in 
order to obtain (17-23). Note, that these equations are regular (for r ^ 0), since 
it holds a 7^ in a neighbourhood of H . This follows from (2) with the identities 
A'^ = a,5{j = e-'^(^) on E. 

In Lemma 3 we shall see that the constraint equations (21) and (22) will propagate 
off from t = . 



Lemma 3. Given as a C^-function and i?H as a C^-function of r G [0,rb]. Let 
;7H,rHeCi([0,rb]) fulfil 

Fh = R'li — Uu 

RhT^ - RhU^ - Rh + 2mH = with 



niH := / 47r [phR'h - (P (Ph) + Ph) ■ Uh] ■ R^ (x) da 
Jo 



(24) 
(25) 

(26) 



on [0,rb]. If (p,R,U,T) is a C^-solution of the evolution equations (17-20) on 
]0,rb] X [-T,T] with 



p = PH, R = Rh, U = Uyl, F = Fh on _ff, 



(27) 
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and if (p, R,U,T) can be extended continuously to r = 0, then (p, R,U,T) also 
fulfils the constramts (21,22) for all t G [-T, T] , r G ]0,rb]. Moreover, R is C'^ 
on ]0,rb] X [— T, T]. Thus (/!>,«'', ds^) of (6), (7) with (11) fulfil Einstein's equations. 

The proof of this lemma is deferred to Appendix A. 

Remark 4- Conversely, if Einstein's equations are fulfilled, then (10), (12), (17-20), (24) 
and (25) are valid for (pn, Rn, Uu, Th) = [{p, R, U, r)]jj. 



4. Constructing black hole spacetimes by solving Einstein's equations 

At first we shall use Einstein's equations (for perfect fiuid) in coordinates, which are 
- in contrast to our (r, t)-coordinates| - regular at the centre of symmetry. We shall 
solve these equations for given data on some auxiliary]: (spacelike) hypersurface F and 
thus obtain a perfect fiuid spacetime Vf with a regular centre. Then we shall construct 
a further perfect fiuid spacetime Vh (as a solution of our (r, t)-equations (17-22)) such 
that [Vf U Vh] is a regular extension of Vf- The spacetime Vh is determined by data on 
a hypersurface H , which will turn out to be a part of an event horizon. We shall finally 
attach to [Vf U Vh] a (unique) vacuum Vq and then gain a 'black hole spacetime'. 



4.1. The perfect fluid spacetime Vf being a neighbourhood of the centreline. 



Lemma 5. Given any smooth 3-manifold F and C°°-data-functions pf,Uf of r ^ [0, ro] 
(some fo > 0). The quantities Py,Uy will become scalar-fields on F (once we have 
chosen suitable coordinates). We assume for pY, Uy 



2n + l 



0, 



f=o 
d?7F' 



dr 



0, for n = 0, 1,2, 



df 



< 



-I f =0 




PY ix) x^Ax I < 1 + (?7f {r) f for r G [0, ?o] 



(28) 
(29) 
(30) 



ion (ds^ , p, u^) of 



Then there exists a unique local (spherically symmetric) C^-solut 
Einstein's equations (for perfect fiuid) on some 4-manifold Dp, where the solution 
assumes the above (essential) data: 

p = Py, U = Uy on F, 

where we use (f , t)-coordinates with ds^ of type (13) fixed hy F = {{r,t) | f = O} 
and i?(f, 0) = r. Moreover, F C -Df is spacelike. Furthermore, there exists a 
neighbourhood Vf C -Df of the centrepomt O, where the auxiliary condition 



{Df + A) {T + U) + ATTpR < on V"f 



(31) 



is fulfilled and where Vf is suitably shaped; i.e. Vf is spherical symmetric and the 
(simply connected) boundary of Vf is spacelike, except for the part between Q and 

t Our (r, t)— equations (17—22) are singular at the centre r = . 

I We shall not give data on the null— cone i7p (see Figure 1), since i7p is a singular surface at the 
centrepoint 0. 
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P (see Figure 1); the points Q and P are joined by an outgoing sound-characteristic, 
whereby P of Figure 1 lies on a null-ray emanating from the centrepoint O ^ F . 

Legend 



line segment from OtoP 

H extension of 

Figure 1. The (spherical symmetric) perfect fluid spacetime Vy U is determined 
by data on F and H . The spacehke hypersurface F contains the centrepoint O. Vji 
is predicted from H \ [neighbourhood of O]. The hypersurface H is null (slope 45°). 

Remark 6. Equation (28) are the necessary and sufficient conditions for C°° -regularity 
at the centre O F . By saying C°° {C^ , respectively) - without referring to 
a particular coordinate system - we mean C°° {C^ , respectively) with respect to 
admissible coordinates. (The polar-coordinates (r,t) are non-admissible at f = 0.) 
Furthermore [R]-p , Uy are C°° with respect to the non-admissihle coordinate f ; but 
they are not (at the centre) in admissible coordinates. However, [-R^]p,[-R]p • 
Uy,<^s'^,Py are C°° in admissible coordinates. As to further details we refer do 
Appendix C. 

Proof of Lemma 5. The data conditions (28-30) are (at t = 0) phrased in (r,t)- 
coordinates. For Einstein's equations however, we use coordinates, which are - in 
contrast to the (f , t)-coordinates (or (r, t)-coordinates) - regular at the centre. The 
equations have a unique (local) spherically symmetric solution (around the centreline) 
for the Cauchy problem; this has been pointed out (by referring to [11]) in that part 
of [12], which deals with the domain containing the centreline (and no part of the star 
boundary). In [12] has been used the actual Cauchy-data (the inner metric of F and 
the extrinsic curvature of F). These data can be calculated (in admissible coordinates) 
from our essential data pY, Uy (see Appendix C). - We obtain (31) from (29) and the 
following lemma concerning the re^a/an^jz-properties of the centreline. 

Lemma 7. Given a spherically symmetric perfect fluid C°°-spacetime V2 - such as Vf 
or Dy (of Lemma 5) - with a regular centreline C . We deflne 

S := {(r,t, ^, ^) I r e [0, ro],t e [-T, T],^=^o,^ = M 
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for some fixed ipo,^o (and sufficiently small ro,T). We consider S as 2-dimensional 
submanifold (of V2), which contains the boundary C at r = 0. Then it holds 

peC\ R,T,U eC°{j:), DrP,Dtp,^ &C^{^), 
R = U = and T = 1 on C; 
for X := (A, drR, DrV) and Y := (DrU, DtU, DrT, DtT) it holds 
lim[X(r,t)],^o = (0,1,0,0) and 

Km [Y(r, t)]^^, = lim [Y(f , i)] , (32) 

where we expressed Y in the right-hand side of (32) in (f , i!)-coordinates used in the 
context of (30). Furthermore, there exists 

lirn [drdrR{r, t)]^-Q ■ 

Proof. We use [6] and the equations (Al, A3), (17-23). 

4.2. The perfect fluid spacetmie Vh hemg a neighbourhood of H \ Nq. 
Let Hp C Vf denote the future null-cone of O G _F (see Figure 1). 

Theorem 8. Given a hypersurfacej H , being any (sufEcently smooth) extension of the 
null-cone Hy (see Figure 1) and given the essential data 

/9h,-Rh on H (33) 
Pb on the star boundary B (pb of (2)) (34) 

where pn, Ra G (H \ {0}) are extensions of [(p, -R)]fjp of spacetime Vf (of Lemma 
5) and where (pn, Rn) fulfil the so-called (pn, -RH)^conditions (47-53), in particular 
the 'black hole' condition (50) (see also Remark 9). From the essential data (ph, Rh) 
the non-essential data can be calculated by (43). Let there be given any (sufficently 
small) neigbourhood #0 of the centrepoint O (on Hp) fulfilling #0 C Vf. Then the 
(r, t)-evolution equations (17-20) have a unique local solution (p, R,U,T) assuming 
the above data 

(/o, i?, [/, F) = (pH, -Rh, t^H, Th) on H \No and p = pb on B. (35) 

The solution (p, R,U,T) defines a perfect fluid spacetime {ds'^,u^,p) with (6, 7, 11) 
in some domain Vh being a neighbourhood of H \ Nq- The scalar fields p,U and 
F are C^(Vh) and R is C^(Vh). Moreover, H (B, respectively) is null (timelike, 
respectively). H will become a piece of an event horizon. Furthermore [Vf U Vh] is a 
perfect fiuid spacetime, which is a -extension of Vf- 

Remark 9. a) The data surface iJ is a null-surface but not a characteristic surface of 
the system of partial differential equations (17-20), since there do not exist spherically 
symmetric gravitational waves. 

b) The data conditions on (pn, Rn) mentioned in Theorem 8 are fulfilled for a 
large class of data (/9h,-Rh) (see Theorem 14). These conditions consist of several 

t see Remark 9 c) 
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mequaltttes on H and several equality relations at the point _Pb (of Figure 1), one of 
these equations at the point _Pb is the Mack hole condition (50). 

c) The evolution equations used in Theorem 8 are written in our (r, t)-coordinates. 
These coordinates have also to be used for equations (33-35); this means H = 
{(r, t) I t = 0, r G [0, rb]} , B = {(r, t) \ r = rb} and equation (35) reads 

Pir,^) = pn{r),...,p{n„t) = pi, for r G [e, rb] , 
where {r,t) = (e, 0) describes the boundary of Nq fl H. 

The proof of Theorem 8 is deferred to Appendix B. 

The construction of the final black hole spacetmie. 

Theorem 10. a) We give - for details see Lemma 5 - data (p-p,U-p) on an auxiliary 
hypersurface F; we gain a unique C°°-solution of Einstein's equations (for perfect 
fluid), where the solution assumes those data p^, IJ-p. The solution is defined in some 
domain Vf (at least); Vf contains a neighbourhood of the centreline. 

b) We give - for details see Theorem 8 - data (pn, -Rh) on a hypersurface H; we 
obtain a unique — solution of the (r, t)-evolution equations (17-20) in some domain 
Vh (see Figure 1), whereby the solution assumes those (essential) data (/9h,-Rh) on 
H being a null-surface, which will become part of an event horizon H . [Vf U Vh] is a 
perfect fiuid spacetime, which is a C"'^-extension of Vf (see Theorem 8). (i?^ is C'^ on 
[T^^U^h]). 

c) We finally attach an (up to extension) unique vacuum spacetime Vq to our 
perfect fiuid spacetime [Vf U Vh] and thus gain a regular and weakly assymptotically 
simple C^-spacetime 

V := [VfUVh\UVo; 

Vo is a Schwarzschild spacetime. We imbed a smooth spacelike, spherically symmetric 
and simply connected hypersurface S (see Figure 2) in V such that S fl [Vf U Vh] lies 
in the past of H and that S ClVo tends to spacelike infinity of Vq - The final spacetime 
V is then a 'black hole spacetime predicted from 5", i.e. 

(i) the smooth (spherically symmetric) null-surface extension H of H (see Figure 
2) is a regular event horizon, i.e. H = J~ (iT""*") ( denotes null infinity, notation 
of [13]); moreover, in the future of H lies a trapped surface area, (which lies in 
parts inside the body,) and a piece of the Schwarzschild singularity.; 

(ii) The future development of S, -D"*" {S), contains a neighbourhood of the event 
horizon H; the closure (with respect to V^ U J') of D'^ (S) contains ^T""*"; 

(iii) D'^ (S) contains a neighbourhood of the centreline and of the star boundary. 

Proof of Theorem 10. The vacuum Vq can be attached at the star boundary B, since 
it holds p = pyj and hence p = on B. Thus the Vq glued together to the perfect fiuid 
becomes a C"'^-spacetime (see [12], which refers to the thesis of Kind). 

Theorem 11. Our spacetime V induces Cauchy data do on S (see Theorem 10). 
From these data it evolves within finite eigentime a unique black hole spacetime 
Vdo := V^ n J'^(S) including an event horizon and trapped surfaces. The above data 
do on S are regular and mnociuous, i.e. there are no points on S, which belong to the 
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horizon or the trapped surface area (see Figure 2). The formation of the horizon (and 
the trapped surfaces) is stable for spherically symmetric variations of the initial data 
do (on S) being in a sufficiently small uniform C^-neighbourhood and for arbitrary 
variations of the equation of state belonging to a C^-neighbourhood. The above data 
do belong to a class of data, which is at least as large as the class of data posed on H 
(in Theorem 10). 

Proof of Theorem 11. We restrict our attention to the matter part of the spacetime, 
since the vacuum is uniquely determined by the matter part. Let us give our data do 
on S (and p = on 5). Then we obviuosly obtain a unique solution [ds^, p] ^ on 
Vdg '■= [Vf U Vh] n J'^(S) assuming the above data. This could be proved basically 
by [12] along with the use of our differentiability classes and coordinates, which are 
analogous to our (r, t)-coordinates. Moreover, this treatment induces stahtUty (with 
respect to spherically symmetric variations): for any data d on S, which lies in a 
sufficiently small uniform C^-neighbourhood of do, the corresponding C^-solution 
[ds^ , m'', p] J lies in a given small uniform C'^-neighbourhood of [ds^, p] ^ . This 
implies in particular the existence of a null-surface Hd with 2md = Rd at Hd H B, 
since Hdg also has this property {Hdg is our previous H , which was the matter part 
of our horizon). The null-extension of Hd (into the attached vacuum) is an event 
horizon. Moreover, the property of trappedness and collapse is also stable. 



5. Admissible initial data 

Our main theorem (Theorem 10) is a special formulation of the following existence and 
uniqueness theorem. Given any set of so-called admissible initial data (pn, -Rh, t^H, Th) 
then the quasilinear hyperbolic initial-boundary problem (17-20) with (35) yields 
a spacetime Vh with the properties stated in Theorem 8 (, provided the data are 
extensions of [{p, R, U, r)]^^). Here we call the functions (pn, Rn, t^H, Th) of r G [0, rj] 
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admissible initial data, if they fulfil 



constraint equations: 

non-negative pressure: 

collapse: 

light escape: 

black hole condition: 

corner condition: 

regular centre: 



Rii {Rii - 2Uh) = 1 - 2mH/i?H, Th = 
(equivalent to (24,25)); 
p(P¥i) > with equality relation in rb; 
?7h < on ]0,rb] ; 

> on [0,rb[; 
2mH(rb) = Ruirh); 
right-hand side of (17) equals at r = 
(for {p, R, U, r) = {pH, Rh, Uh, Th) ) ; 
i?H(0) = 0,i?ij(0) = l,[/H(0) = 0. 



Rii - Uh (36) 



(37) 
(38) 
(39) 
(40) 
(41) 



(42) 



In general, one can prescribe any two of the four functions (pn, Rh, Uh, Th) and 
obtain the others from the constraint equations. In practice, the most simplej way 
seems to be to prescribe pn and Rh, which we therefore called essential data in the 
previous sections. The non-essential data Uh and Fh can be calculated from the 
constraints as 



Uh=U {ph, Rh) and Th 
where the operator U is defined by 

U{ph,Rh) ■■ [0,oo[^M, r ^ 
with qii := R 

X{r) := 



R'h ~ Uh, 



(R'nY 



Rnir) JO 
H -r iirpHRn and 
qh + 4:Ttp(ph) Rh 



?H • Rh 



R'tt 



dx. 



(43) 

(44) 
(45) 
(46) 



Since we shall prescribe the essential data (pn, Rh) and calculate the other data 
by (43), we shall now rewrite the conditions (36-42) as conditions on pn and Rh- 

Lemma 12. Let pn G and Rh G for r £ [0,rb]. Then {pn, Rh,Uh,^h) 
with Uh G C'"'^([0) ''b]) are admissible initial data (i.e. fulfil (36-42)) if and only if 
Uh = U {ph, Rh) , Th = R'h ~ Uh and the following {pn, Rh) — conditions are fulfilled: 



p> Ph on [0, rb[ p{rh) = Ph 
U{ph,Rh) < in ]0,rb] 
i?H > on [0,rb[ 
i?ij(rb) = 

9h (j'b) = Stt ■ pb (ph, Rh) (r-h) 
qnirh) = 

i?H(0) = 0,i?ij(0) = l 



(non-negative pressure), 

(collapse), 

(light escape), 

(black hole condition), 

(corner condition), 

(see Remark 13), 

(regular centre of symmetry). 



(47) 
(48) 
(49) 
(50) 
(51) 
(52) 
(53) 



Before we proof this lemma, let us exhibit a few properties of these {pu,Rb)~ 
conditions. 

t For given (phj^h) the quantity [C/jj * ^h] is a solution of a linear differential equation being 
equivalent to the first constraint equation in (36). 
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Remark 13. a) Note, that (48,49) are inequalities on H , whereas (50,51) are equality 
relations at the point H Cl B. The corner condition (51), which is equivalent to (41) 
is a technical condition ensuring a convenient degree of differentiability. 

b) The inequalities (47-49) are stable under small variations (of (ph, Rb)) keeping 
the black hole condition (50) and p(rYj) = pb fixed. This follows from the continuity 
of the operator U . 

c) The black hole condition (50) is equivalent with [^J^^^ = 1, provided it holds 
(36, 48,49). 

d) Equation (49) means (Df + Dt)R > Q on H \ B, where H is - among other 
things - a future-lightray. In general, it holds for any connected future-lightray 
L C [Vf U Vh] the following: L escapes^ to ^T""*", tf and only if (Df + Dt)R > 
on L \ B, provided L hits the star boundary B and p > on Vf U Vh U 5, 
(which holds for sufficiently small Vf and Vh.) The escape-condition can be proven 
by [6] and the fact that {Df + Dt)R is strictly positive (negative, respectively) on 
Bf^Vof^[J-{H)\H] {Bf^Vof^[J+{H)\H], respectively) . 

e) Equation (49) is equivalent with [^]^^^ < 1, provided it holds (36,48,53). 
Furthermore, (49) prevents F = 0, which would make the evolution equations (17-20) 
singular. Moreover, F < would mean a 'neck' (see [5]). 

f) Equation (52) implies that the integrand of x of (46) is - despite (50) - 
continuous. Thus UY{{p]i, Rn) is also at r = rb (cf (46)). Conversely, the first 
constraint equation of (36) and i?^(rb) = 0, which is a consequence of (36-42), imply 
(52). 

Proof of Lemma 12. At first we show the equivalence of the constraint equations (36) 
with (43) by rewriting the first equation of (36) as 

R'u RnUn = mn + Rn- ({R'^f - l) /2 (54) 

with = 47r(pHi?H- (p(PH) + /0H);7H)i?H of (26). Then [^(54)] / R!^ becomes 
a linear differential equation in the new unknown W := [i?H • t^n]- For initial value 
11/(0) = it has a unique solution W , which is explicitly writable as an integral. 
Indeed, we defined in (44) our U {pn, Rb) as W/Rn; this yields the 

equivalence of (36) with (43) on ]0, rb[ at least. (55) 

If (36-42) are fulfilled with Un e €'^[0, rj) then Uh(0) = = U (ph, Rh) (0) by 
Hospital's rule. Because of = (U {pa, Ra))' on ]0, rb[ there exists 

lim (p-ii,Rb))' (r) = lim 

Together with lim^^rb (^ (pn, Rh)) (f) = Unlryj) g]0, cxd[ it follows the existence of (a 
finite) lim^^rb x'(^) = x'(^b)- This implies qii(rYj) = due to (50), which follows from 
(36), (40) and (R'^ - 2Uh) (rb) ^ 0. 

Conversely, if it holds (43) and (47-53) then it follows 

?7h(0) = 0, (36) for r = and Un G ([0, rj) 

from (44), where we use Hospital's rule at r = (r = rb, respectively) bearing in 
mind (53) ((52), respectively). From this and (55) we gain (36) on [0,rb[, which can 

+ If L hits B n then H is the unique extension of L, i.e. H 'escapes to 3^'^^ in the sense that 
there exists an affine parameter, which tends to infinity (when moving towards the future). 



i?H 



U (pH, -Rh) 



2Rh 



(r). 
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be extended to [0, rb] by the continuity of the quantities in (36). Then the black hole 
condition (40) follows from (50). 

It remains to show the equivalence of the different formulations of the corner 
condition. The evaluation of (41) just yields (51), where we used (17) for (p, R, U, V) = 
{pn,Rn,Un,Yn), a = A'^,B'^e^^P^ = l,p= 0, (22), (40), and (43,46)' with (50) and 
the rule of I'Hospital for Qh/R'h- This finishes the proof. 

We shall now show that the (pn, -RH)^conditions are not too restrictive. At first, 
we choose pn as free tmttal data (with (47)). Then the amount of freedom left to the 
choice of i?H amounts to the cardinality of ([0, 1] ,]0, — oo[) (at least). 

Theorem I4. Let pn be given with (47) and 
167r I max{pii(r)} I < 1. 

\[0>'-b] / 

Then for any auxiliary C^-function h : [0, 1] ]0, — cxd[ we can construct an 
'admissible' i?H = -Rh such that (/9h,-Rh) fulfil the (pn, -RH)^conditions (47-53), 
i.e. (ph, -Rh, t^H, Th) are admissible initial data, where (J/hiTh) are defined by (43). 
Furthermore, our construction yields R'^ ^ R}^ for h ^ h. 

We shall describe the main ideas of our proof. Since we are concerned with 
quantities (such as pn, -Rh, t^H, • • •) defined on H only, we shall drop the index H 
and bear in mind that these quantities depend on r only. From (44) one sees that 

q = R" + A-KpR < , (56) 

R' <1, (53) and (49) are sufficient for [/ < 0. 

Therefore we do not choose R with (48), (49) etc., but we choose a negative q Q C° , 
which defines R := TZ (q) as the solution of 

R" + iTrpR = q with i? (0) = and i?' (0) = 1. (57) 

The operater TZ : C° ^ is continuous if one endowes the range of TZ by the 
C^-norm (uniform topology). 

In order to gain an R := TZ(q) as in Theorem 8 we choose g to be g = go on a 
sufficiently small part of Hp, i.e. for r G [0, e]; qo is defined as 

qo := [left-hand side of (31)]^^ < 0. (58) 

Moreover, we postulate q to be near rb, since (/'(rb) occurs in the corner condition 
(51). 

Lemma 15. Let p and h be given as in Theorem 14. Then a real parameter a and a 
positive parameter r exist such that 

9 := 9(<7,r) := Kr [qh + T]a) (59) 
induces (in the sense of Theorem 14) an admissible R := TZ (q) (TZ is defined around 
(57)). Here we define qh as a continuous function with qh(r) := h ^ rV-^Se ) 
r e [3£,rb - 3£], qu = qo on [0,£], < on [0,rb - e[ and = on [rb - £,oo[, 
where qo was defined in (58). The auxiliary functions and rja are defined below. 
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Remark 16. The function of Lemma 15, which determines R := TZ(qh), carries the 
information of h. Moreover, if the assumptions of Theorem 8 are fulfilled, then R := 
T^{<lh) becomes an extension of by choosing go := [left-hand side of (31)]^^ < 0, 

go G C'^([0,£]). Furthermore, is near r\^. 



Scetch of the proof of Lemma 15 including the definitions of Kj and rja. 

The parameter r (cf (59)) is used for positioning the (first) 
zero of R' = \TI {q{a,T))\ at r = r\^ 



(60) 



When we increase r (decrease, respectively) then the zero of R' decreases (increases, 
respectively). Kj is defined as = 1 + (r — 1)k, where k : [0,cxd[ [0,1] is any 
fixed C°°-function with support C [s, rb — s] and k = 1 on [2e, rb — 2e]. We define 
a zero-searcher ( such that ((q) is the first (positive) zero of the derivative o{TZ(q), 
provided q lies in a (sufficiently small) neighbourhood of a negative function! . 
It can be shown that ri , T2 exist with 

C (ktjS) < J-b < C (ktiS) and n < t-j (61) 

for all g of a C'^-neighbourhood of qu- Moreover, the operator ( is continuous there. 
Let such a g be given. Then we can find (according to the intermediate value theorem) 
a Kr (with r G [ri, T2]), which does its job: C('^r • ?) = ''b- 

The other family of auxiliary functions rja controls the derivative of q(a,T) at 
rb with the aim to fulfil the corner condition (51). Unfortunately, a variation of 
q' (rb) = rj'^ (rb) = cr violates the requirement g < and also handicaps the aim 
concerning the choice of r (see (60)). The first item does not matter provided \ria \ is 
small enough, say \ria \ < 6q for some constant 6q < e. 

We define ri„ : [0, cxd[ [—^q , ^q] by 

T]a (r) := 6a sign((T) • fj — \a\ 



with := min {6q ,s\a\}, where fj : [0, cxd[ ^ [— 1, 1] is any C"'^-function with fj' (0) = 1 
and support C [—1,1]. 

It can be shown that {</ G \ <! = Ih + i1<j with cr G M} lies in a neighbourhood of 
qh such that there exist ri, T2 with (61). We define the set of 'black hole'- parameters 

E ■.= {(a,T)emx [n, r2] I i?'(rb) = for i? := 7^ (?(,,,)) (cf (59))} 

(cf condition (50)). For any given cr G M the set E contains one (c, t) at least; 
thus (a, t) ^ E means no restriction on a. Up to now, we know that (p, R) with 
R := TZ{q(a,T)) fulfil the (pn, -RH)^conditions, if (c, t) G E, provided the corner 
condition (51) is satisfied for qn = q(a,T) and Un = U(p, R). Due to t)(^^->) ~ 
corner condition is fulfilled, if it holds 

a = 87rpbUhi(T,T), (62) 

where the function J/b : -E" ^ M is defined as [/^{a, r) := [U {p,TZ (?((7,r)))] (''b)- Thus 
it remains to show the existence of such a (a, r) G E, which fulfills (62). Since Uyj is 



t Then {7Z{q))' is strictly monotonically decreasing, which makes C well-defined and continuous. 
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uniformly bounded and continuous]:, we have to look for these (c, t) G E with (62) 
only in some compactum [(Ti,(T2] x [ri,r2]. We show that such a (c, t) exists by the 
intermediate value theorem along with some technical arguments. Thus we have found 
an admissible R = TZ {q(a,T))- 

It is not hard to see that ^ qj^ lead to admissible R^ ^ R^ . 
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Appendix A. 

Proof of Lemma 3. For r 7^ it holds that a C^-solution (p, R, U, F) of (17-20) is also 
a C^-solution of 

dtp =-e^W(p + p) ^^"'^+"^ 

dtR = e^^P^U 

A'^ dp 
(to be proved via (15) assuming (2)). 

Remark. Note that (Al) turns out to be Misner's evolution equations (with ^ 
chosen as j4^), provided one formally substitutes his — 2m/i? by V'^ — U'^ — 1. 

We now shall need dfdfR = dfdfR and hence prove i? G by defining a scalar 
field R (unique up to an integration constant) with 

''r-R = 1 P2 • r and dtR = e^(^) U , (A2) 

Att ■ n(p) ■ R^ 

(i^T^Pp?^ ' ^) = (e^^'P'>U^ because of (Al). Thus (A2) implies R e . 
From (15) we get for r ^ 

^^^=i-f^-r + i^He^(^)^, (A3) 

thus ^drR^ ^ = i?^ for r 7^ because of (4,5,24). This and fixing R by R{ro, 5H(ro)) = 
Rii(ro) (for any fixed rg G ]0, ry^]) gives R = R on H (for r ^ 0), which (together with 

+ In order to proof this we define tlie continuous operator 7" : | / G near | /(^'b) = o| ^ 
as (r(/)) (r) := Then the continuity of S ^ C«, (a,r) ^ x' = ^^''"r^Tz'^r^))?/'"^'^ 

can be shown. 
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dtR = dtk) yields i? = _R e in ]0, r^] x [-T, T]. From (A3) and i? = _R it follows 
(21). 

Analogously to the construction of R, we define fh Q (unique up to a constant) 

by 

A' 

drm = Vr ■.= 4:TTpR^ -^^ - T and dfrh = Vf := -'iirpR^Ue'^^P'' (A4) 

47r ■ n(p) ■ R^ 

for r g]0, Tb], (t — Bii(r)) G [—T, T], where dfVf = dfVf, which follows from dfdfR = 
dfdfR and (Al). Now we consider the scalar field rh in (r, t)-coordinates; we use 
(15),(A4),(21) and obtain for r ^ 0: 

drfh = drin, (m of (23)) (A5) 
m(r, 0) = m (ro, 0) + / drm(x, 0)dx, (A6) 



m(r,t) = m(r,0) + K*"- ^^](r,t)=(r,:.+i3H(r)) (A7) 

Furthermore, drm(r, 0) = is continuous on [0, rj (express again Sr-R by F of (21)); 
this and (A6) implies that limr^o rh(r, 0) exists; now, rh was fixed up to a constant; 
we choose this constant such that limr^o rh(r, 0) = 0. This and the fact that Vf tends 
continuously towards 0, implies limr^o rh(r,t) = (using (A7)). From this and from 
limr^o m(r,t) = (m of (23)) we gain the identity m = m by (A5). 



Equations (Al) and (A4) along with rh = m implies dj [F^ — [/^ — 1 



2m 1 
R I 



-e'^(^)g • [T^ - U^ -1 + ^]. This and [r^ - [/^ _ i + 2m j ^ ^ q ^^^^ ^25)) 
implies that the constraint (22) is valid. 

Due to (A4) Misner's mass (cf [5]-[7]) can be identified as fh, i.e. our m and 
Misner's m are the same scalar fields. This, (Al) and the constraints (21,22) imply 
Misner's equations (with ^ chosen as A'^), as stated in the remark under (Al). Thus 
Einstein's equations hold. - Note, that our evolution system (Al) appears in a form, 
to which Wesson and de Ponce do not object in [14]. 

Appendix B. 

Proof of Theorem 8. Courant and Hilbert [15, p 425, p 476, pp 461 ffj proved 
the existence and uniqueness of the solution of the following mixed problem for a 
quasilinear system (in two variables). Let u = (R, p,U,T) be the vector of the 
unknowns. Then one can write (17-20) as 

dtu + M- dru = f 

where the 4 x 4-Matrix M and / depend on (r,t,u(r,t)) . The eigenvalues of M are 
Xi 2 = and 

±,/^47rni?2 e'^(^) 
, V dp 



We can choose as associated left eigenvectors 

= (0, 1, 0, 0),V2 = (0, 0, -U, F) and ^^3,4 = IJ^' 0' ±^7"^' 
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i.e. XiVi = Vj- • M for i = 1,2,3,4. Our system is hyperbolic, since these eigenvectors 
are linear independent for u = mh- The canonical functions gi := vi • u = R and 
g2 := V2 • u = V'^ — propagate along the stream lines, while 



propagate along out-/in-going sound characteristics (cf [16]). So our mixed problem 
is well-posed [15, p 472] if we give initial data u (r, 0) = mh (f) and boundary values 
for §4, namely 



These data mean p{r\-j,t) = pb (for t G [~T,T]), i.e. vanishing pressure on the 
boundary of the star. This initial-boundary value problem has a solution with 
Lipschitz continuous first derivatives, provided M, /, A,-, Vi (and mh respectively) have 
Lipschitz continuous first derivatives with respect to t,r,u (and r respectively). 

We have to show the uniqueness of the solution in Vf H Vh, i-e. that part of the 
spacetime in figure 1, which is marked twice (by parallel lines in horizontal as well as in 
vertical direction). The boundary of Vf H Vh is such that the slopes of its tangents are 
not steeper than those of the sound characteristic, which is the fastest characteristic. 
This yields 



where denotes the (future and past) domain of dependence with respect to the 
sound characteristic. On the other hand, we have solutions of the evolution equations 
(e.g. (Al) for r > 0) together with the constraint equations (21,22), since Lemma 3 
is extendable to any simple connected domain with simple connented stream lines by 
extending m as m. These solutions of Einstein's equations exist in Vf as well as in 
Vh and coincide on iJ fl Vf H Vh. This and (Bl) yields coincidence of the solution in 
Vf n Vfj. This finishes the proof of Theorem 8. 

Appendix C. 

Scetch of the proof of Remark 6. It has been pointed out in [12] that (28) is a necessary 
condition for C°° -regularity at the centre. We show that these conditions are sufficient: 
the actual data of our Cauchy-problem in Lemma 5 are the inner metric of F 
(a, j3 runs 1,2,3) and the extrinsic curvature of F . They can be calculated from 
our (essential) data pY, Uy as 






[V^fH^h] C D,{Hi), 



Hi:=Hn [Vf n V^h] 



(Bl) 



9al3 



E - 1 




(CI) 



with 




(C2) 



where we used (13), R{r, 0) = f and the coordinates! 

xi = rcosLpsm{), X2 = rsunpsmi), xs = rcosi9 

t These coordinates turn out to be admissible, provided it holds (28) (see also [12]). 
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in order to obtain (CI) (as to (C2) see [12]). 

We turn our attention to the regularity at the centre f = 0. Due to (28) 
we can extend (Uf, repectively) as an even (odd, respectively) (7°°-function of 
r G [—To, To]. This and (C2) implies that E can be extended as an even C°°([— fo, tq])- 
function. Thus there exists Ei G C°°([0,ro]) with _E'i(f^) = E{r). From this and 
= E{{)) - 1 = £;i(0) - 1 it follows the existence of some E2 G C°^([0,r^]) with 
{Ei{P) - 1) . = E2{P) = E2{xl ^xl^ xl). This, E{r) = Ei{P) and (CI) imply 
^ ^® functions of x^ and hence as a tensor field on _F, which also holds for 
A'^^ and (lo be shown analogously). 

This implies that the Cauchy development of the (7°°-data (^^^, A'^^, pp) is 
(see [12], who used [11]). 
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